Binding of a mobile hole by an impurity potential in the t-J model: parity breaking 
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We revisit the problem of a single hole moving in the background of the two dimensional Heisen- 
berg antiferromagnet. The hole is loosely bound by an impurity potential. We show that the bound 
state is generically a parity doublet: there are parametrically close bound states of opposite parity. 
Due to the degeneracy the bound state readily breaks local symmetries of the square lattice and 
this leads to formation of the long range spiral distortion of the antiferromagnetic background. A 
direct analogy with van der Waals forces in atomic physics is discussed. 

PACS numbers: 74.72.-h 75.25.+Z 75.30.Fv 



I. INTRODUCTION 

The problem of single hole binding by an attractive 
potential in the t-J model is of fundamental importance. 
In the physics of doped Mott insulatoris and in particular 
in the physics of the cuprate superconductors, the system 
play a role analogous to the that played by the hydrogen 
atom in atomic physics. 

We have in mind, for example, La2Cu04 with a La 
ion replaced by Sr. Alternatively it may be Ca2CuC>2Cl2 
with a Ca ion replaced by Na. An important point is 
that the attractive center (Sr ion in La2Cu04 or Na ion 
in Ca2CuC>2Cl2) sits in the center of a square of four Cu 
sites. This means that the attractive potential itself does 
not break the local square lattice symmetry. There are 
various aspects of the bound state problem: the symme- 
try/parity of the bound state, the structure of the spin 
fabric, and in the end the particular value of the binding 
energy. We argue that the symmetry issues are the most 
important ones. 

There have been several studies of the bound state 
problem. These are mainly small cluster exact dioganal- 
izationsi^i^. A generic limitation of this approach is 
the small cluster size and as a consequence sensitivity to 
boundary conditions. In spite of this limitation a very im- 
portant observation has been made already in the early 
work^: the ground state is almost degenerate with an- 
other state that has opposite parity. Dependent on the 
parameters of the model, the ground state belongs ei- 
ther to the two dimensional E-rcprcsentation of the C^ v 
symmetry group of the Hamiltonian or to the A\ rep- 
resentation^. However, there is always a very low-lying 
excitation of opposite parity. 

A semiclassical solution of the bound state problem 
was obtained in Ref£ Generally a semiclassical approach 
can be justified in the limit of a large radius of the bound 
state. According to the semiclassical solution the bound 
state generates a long range (oc 1/r) spiral distortion of 
the spin fabric as it is shown in FigfT] The figure shows 
staggered spins. It is obvious from Fig. Q] that the semi- 
classical solution does not have a definite parity. It does 
not fit in any representation of the symmetry group of 
the square lattice. So the solution spontaneously breaks 
the local square lattice symmetry. On the other hand 
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FIG. 1: (Color online). Distortion of the staggered spin fab- 
ric (small arrows) by the Sr-hole bound state. The left pic- 
ture corresponds to the pseudospin directed out of the page 
and the right picture corresponds to the pseudospin directed 
in the page. Shaded area corresponds to the hole localiza- 
tion region. At large distances spins are directed along the 
orthorhombic b-axis due to pinning by Dzyaloshinksi-Moriya 
and XY anisotropics. 



it is impossible to have a spontaneous violation of an 
exact symmetry of the Hamiltonian in a finite system. 
Therefore, the exact parameter that justifies the semi- 
classical solution^ has remained unclear. The purpose of 
the present work is to elucidate this parameter and hence 
to elucidate the physical meaning of the solution with vi- 
olation of exact symmetries. We show that the physics of 
this system is similar to the physics of a hydrogen atom 
in an external electric field. The parameter that justifies 
the semiclassical solution is the small splitting between 
states of opposite parity. The splitting scales as the bind- 
ing energy squared and hence it is infinitesimally small 
for a shallow bound state. 

The structure of the paper is as follows. In Section 
II we explain the analogy with the hydrogen atom. In 
particular, we consider the conditions when a long range 
tail of the dipole electric field can be generated by the 
atom. In Section III we briefly review known properties 
of an unbound single hole moving in the antiferromag- 
netic background of the t—J model. Section IV addresses 
the limiting case of very strong binding. Here we present 
results of exact diagonalizations for the 4x4 cluster em- 
bedded in an antiferromagnetic background. In Sections 
V and VI we consider the weak binding limit and discuss 
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symmetry properties of the bound states. Section VII 
addresses the parity breaking and generation of the local 
spin spiral. In Section VIII we exclude a possibility of 
the local Charge Density Wave (CDW) formation. Our 
conclusions are presented in Section IX. 



II. HYDROGEN ATOM 

Consider a hydrogen atom in the ground Is state. The 
size is about one Bohr radius, ag. Since the atom is 
neutral the electric field at distances r>ag decays ex- 
ponentially. Let us consider the same atom in the n = 2 
state, either the positive parity 2s-state or negative parity 
2p-state. Importantly, they are degenerate. Because of 
the degeneracy an infmitcsimally small external electric 
field, E ext — > 0, will mix the opposite parity states 



1> 



1 



V2 



2s) 



1 



V2 



2po> 



(1) 



Here \2po) is the state with zero projection of the angular 
momentum in the direction of the external electric field. 
The state ([T]) possesses a static electric dipole moment 
d ~ eas- Hence a static dipole electric potential and a 
static electric field are induced outside the atom, r 3> as 

(d-r) „ , , d 3(d • r)r 

¥w(r) = , E ind (r) = -- + -L-gJ- . (2) 

Due to the small but nonzero energy splitting A between 
2s- and 2p-statcs (Lamb shift), one needs to apply a small 
but finite external field, E ext > A/d, to create the mixed 
state (JXJ) and hence to induce the dipole field (|2|) . Impor- 
tantly, the induced field @ is much larger than E ext . 

One can also look at the problem from another point 
of view. Consider two hydrogen atoms each in the n = 2 
state. The attractive potential between the atoms has 
two distinct regimes depending on the distance r between 
the atoms. The characteristic distance is defined by 
the condition d 2 /r\ ~ A. If <C r <C ^a the potential 
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In this regime the electric dipole fields of the two atoms 
lock to each other. At r 3> the interaction scales 
as 1/r 6 , this is the usual van der Waals regime that is 
due to fluctuating dipoles. We will argue below that the 
dipole distortion of the spin fabric in the two-dimensional 
(2D) t-J model shown in Fig. [1] is fully analogous to Eqs. 
dU),©. The power in the 2D case is different, — ► -5. 
A more important difference is that in the t-J model the 
ground state itself is a parity doublet. The splitting in the 
doublet is paramctrically small at small binding energy, 
it scales as A oc e 2 , where e is the binding energy. 



III. A FREE HOLE PROPAGATION IN THE 
t J MODEL 

The 2D t — J model was suggested two decades ago 
to describe the essential low-energy physics of high-T c 
cuprates£i£i£. In its extended version, this model includes 
additional hopping matrix elements t' and t" to 2nd and 
3rd- nearest Cu neighbors. The Hamiltonian of the t — 
t 1 — t" — J model on the square Cu lattice has the form: 

H t -j = -* C ^' CT ~ ° ] ™ C i' a ~ l " /Z Aa c i"° 

(ij)a ('»j')< 7 (ij")<r 

Here, c\ a is the creation operator for an electron with spin 
a (a =T,I) at site i of the square lattice, (ij) indicates 
1st-, (ij 1 ) 2nd-, and (ij") 3rd-nearest neighbor sites. The 
spin operator is S, = \^ ia o a pc i(3 , and Ni = Y^aAa^a is 
the number density operator. In addition to the Hamil- 
tonian ((T|) there is the constraint of no double occupancy, 
which accounts for strong electron correlations. 

The values of the parameters of the Hamiltonian (|4|) 
for cupratcs are known from neutron scattering, Raman 
spectroscopy, and ab-initio calculations. For La2Cu04 
the values ar e 9 ' 10 i n : 



(4) 



J 
t 
t' 
t" 



140 meV -> 1 , 
450 meV , 
-70meV , 



35 meV 



(5) 



Hereafter we set J = 1, hence we measure energies in 
units of J. In the present work we study generic proper- 
ties of the extended t — J model. Therefore we will vary 
parameters t, t' , and t" in a broad range. 

At zero doping (no holes), the t-J model is equivalent 
to the Heisenberg model and describes the Mott insulator 
La2Cu04. The removal of a single electron from this 
Mott insulator, or in other words the injection of a hole, 
allows the charge carrier to propagate. 

The properties of a free single hole in the t-J model 
are very well studied numerically: see Ref.— for a re- 
view. At values of parameters corresponding to the 
cuprates the dispersion of the hole dressed by magnetic 
quantum fluctuations has minima at the 'nodal points' 
q = (±7r/2, ±7r/2) see Fig. [5J The typical value of the 
quasiparticle residue at these points is Z sa 0.3. By 
changing the sign of t' and t" one can shift the dis- 
persion minima to the 'antinodal points' (±7r, 0) and 
(±7r, 0). This situation corresponds to the electron doped 
cuprates. We will argue below that the properties of 
the shallow bound state are most interesting and rich in 
the regime when the minima of the dispersion are at the 
nodal points. This is the regime which we consider in the 
present work. The dispersion of the hole dressed by mag- 
netic quantum fluctuations is quadratic in the vicinity of 
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FIG. 2: Magnetic Brillouin Zone with a- and b-minima of the 
hole dispersion 



qo, 



e (P) « -ftp? + 0&P2 , 



(6) 



where p = q q . We set the lattice spacing to unity, 
3.81 A — > 1. In Eq.© pi is directed along the nodal di- 
rection and P2 is directed along the face of the Magnetic 
Brillouin Zone (MBZ), see Fig. [2] At values of the hop- 
ping parameters presented in Eq. ([5]) the inverse masses 
are!2, ft « ft « 2.5. 

It is instructive to consider also the weak coupling limit 
( C J, f' = t" = 0. In this limit the quasiparticle residue 
is close to unity, Z = 1 — 0(t 2 /J 2 ), while the dispersion 
reads^i 

e q = 4t^/ j, cos g x cos g y + 2t" f f (cos 2g x + cos 2q y ) , 
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(7) 



IV. 



HOLE BINDING IN THE STRONG 
COUPLING LIMIT 



We include a site-dependant potential attraction to an 
impurity, 



H — H t -j + Hjj , 
H v = = U l c\ <J c il 



(8) 



A very important point is that the potential Ui is sym- 
metric around the center of a plaquette, see Fig. [3j Note 
that the Hamiltonian © is written in terms of electrons. 
Repulsion for electrons [7, > corresponds to hole at- 
traction. Concerning the potential we will consider two 
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FIG. 3: (Color online). The strong binding limit. Attraction 
to the potential center ( a filled red circle in the center of the 
plaquette) is so strong that the hole (small empty circle) can 
hop only within the few sites around the attractive center. In 
the exact cluster diagonalization we assume that the perimeter 
spins (blue arrows outside of the dashed square) are static with 
(S z ) = ±0.3. 



possibilities. The first possibility is the short range po- 
tential which is nonzero only at four nearest sites, see 
Fig.H 



Ui 



= uj2k 



(9) 



where j runs over four nearest sites. The second possi- 
bility is the long range Coulomb interaction 



Q 



(10) 



where Q is a dimensionlcss charge which sits at a distance 
one lattice spacing above the plane. 

Let us consider first the local potential (J9j) in the strong 
coupling limit, U 3> t. The solution of the bound state 
problem in this limit is qualitatively clear. There are 
degenerate states with S z = At each value of S z 

there are states of positive and negative parity. Let us 
consider the lowest bound state in each parity sector. 
Note that there is only one lowest state in each sector. 
This is contrary to the common wisdom that the negative 
parity states are doubly degenerate due to the symme- 
try of the square lattice (E-representation of C^). The 
point is that we consider hole binding on the spin back- 
ground with spontaneously broken SU(2) symmetry. In 
combination with the impurity potential this breaks the 
symmetry of the square lattice and hence destroys de- 
generacy of the negative parity states. Binding energies 
of the lowest bound states are 



-4 



(11) 



Here the sign ± denotes the parity of the bound state. 
We define the binding energy e in the standard way: this 
is the energy of the bound state taken with respect to the 
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minimum energy of a free hole. So e is always negative. 
To find the parity splitting A in the strong coupling limit 
we have performed exact diagonalizations of the t — t' — J 
model on the 16 site cluster shown in Fig. [3] We have 
already pointed out that it is qualitatively important to 
perform the diagonalization on the state with sponta- 
neously broken SU(2) symmetry. Therefore we put the 
cluster in the environment of static perimeter spins shown 
in Fig. [3] by blue arrows outside of the dashed square. 
The magnetization of each static spin has the Heiscn- 
berg model value, (S z ) = ±0.3. Values of the splitting 
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TABLE I: Exact diagonalization of the 16-site cluster (Fig. 
The ground state parity doublet energy splitting and the rms 
charge radius of the ground state for several values of t and 
tl and for two values of the confining potential U. According 
to Eg. A > corresponds to the negative parity of the 
ground state and A < corresponds to the positive parity of 
the ground state. The bound state results for values of t and 
t' that correspond to the free hole dispersion minima at the 
nodal points, (±7r/2, ±7r/2), are presented by large font. 

A within the parity doublet obtained by the 16 site clus- 
ter exact diagonaization are presented in Table I. In the 
same table we present values of the rms charge radius 
of the lowest bound state. The results are presented for 
several values of t and t'. The value of t" in this calcula- 
tion is zero: the cluster is too small to account for long 
range hopping. We have performed the calculation for 
two values of the confining potential, U=0 and U=10. In 
an infinite system the case U=0 certainly does not cor- 
respond to any binding. However, for the cluster, due 
to the imposed boundary conditions, the case U=0 de- 
scribes a well localized state of the hole; in this sense it 
is bound. Our numerical results qualitatively agree with 
those of previous publications^^. A detailed quantita- 



tive comparison is not possible because the previous pub- 
lications have considered spin symmetric clusters while 
we impose a spontaneous violation of the SU(2) symme- 
try via boundary conditions. 

The dispersion of a free hole for various values of t and 
t' is well known from previous wor k 12 i 13 ' 14 . The bound 
state results in Table I for values of t and t' that cor- 
respond to the free hole dispersion minima at the nodal 
points, (±7r/2,±7r/2), are presented by large font. For 
all other values of t and t' the dispersion minima are at 
the antinodal points, (±7r,0), (0,±7r), or at the T-point, 
k = 0. The corresponding bound state results in Table I 
are shown by standard font. Results presented in Table 
I lead to the following observations. 

1) Values of the parity splittings for the 'nodal cases' 
(large font) are very small compared to typical scales in 
the problem. This conclusion is in agreement with pre- 
vious observation^. 

2) Values of the splittings for other cases ('antinodal' and 
T-point') are substantially larger. 

3) For each particular 'nodal' set of t and if the splitting 
A for U=0 is systematically smaller than that for U=10. 

The first and second observations indicate that the 
bound state is close to the parity degeneracy in the case 
of the nodal minima of the free hole dispersion. This 
is why in the present work we concentrate on this case. 
According to the third observation the parity splitting 
is rapidly decreasing when the radius of the bound state 
increases. In following section we consider shallow bound 
states of large radius to confirm these conclusions. 



V. HOLE BINDING IN THE WEAK COUPLING 
LIMIT: THE LEADING APPROXIMATION 

Let us look at the binding problem in the weak cou- 
pling limit, U — > 0, e — > 0. In this case the shallow bound 
state can be built with a hole either from the a- or the 
b-valley of the dispersion, see Fig. [2l Hence the bound 
state has the valley index and the corresponding wave 
function reads 

Vw(r) = e^x^rv) , (12) 

where a = a,b shows the valley and a =|, | shows the 
magnetic sublattice along which the hole is propagating, 
r a is position on this sublattice. Note that here we have 
in mind a real propagation. There are also virtual hop- 
pings of the hole to the opposite sublattice. These vir- 
tual hoppings lead to formation of the free hole disper- 
sion that was discussed in Section III. The z-projection 
of the hole spin is S z = — a . The oscillating exponen- 
tial dependence e l i<*' r ° in (Tj"2"]) is due to the momentum 
q a = (±7r/2, ±7r/2) that corresponds to the valley mini- 
mum. The very smooth function x(r) exponentially de- 
caying at infinity is due to the hole binding to the poten- 
tial. In the case of the Coulomb field, Eq. pI7|) . the wave 
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function is^ 



l Kp -^\e\(rl/0 1+ ri/0 2 ) 




(13) 



The components r\ and T2 in Eq. (|13p are projections of 
r on directions 1 and 2 corresponding to the particular 
valley, see Fig. [5J In the case of the local attraction 
Eq.©, the wave function is 



X 



rllfo)), 



(14) 



where Kq is the Bessel function of the second kind. Note 
that for both Eqs. (fT3|) and (fT4|) the root mean square 
radius of the bound state scales as 



1 



(15) 



It is easy to see that a change of sign of q a leads only 
to a common phase factor in the wave function (|12[) . so 
this is the same wave function. There are only two dis- 
tinct possibilities: the a-minimum, q a = (tt/2,tt/2), and 
the b-minimum, qb = (w/2, —w/2). Thus, there are two 
degenerate quantum states for each value of S z - 

We put the potential center at the origin of the coor- 
dinate system. Then, according to Eq. (fl~2"| and Fig. [3] 
(in this case one has to remove the cluster boundary and 
extend the figure up to infinity) the | f) wave functions 
read 



0T<» 
0? ( r ) 



x 1 



VI = o + n 



- m , 



(16) 



where both m + n and m — n are integer and even. Sim- 
ilarly the | I) wave functions are 

ifw Xb(r) , 



x l 



(17) 



We put an additional factor —i in rp a ^ and ip^i to make 
these wave functions real. Under the parity operation 
x — > — x and y — » —y, the function x( r ) does not change. 
Therefore parities of states (fit))) and (fTT|) are determined 
by the phase factors, and the parities are 
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Thus, in the leading weak coupling limit approximation, 
e — * 0, the ground state is a degenerate parity doublet for 
each value of S z . This explains why values of the parity 
splitting presented in Table I by large font are very 
small. The next Section is addressed to the mechanism 
that lifts the exact parity degeneracy. 



VI. HOLE BINDING IN THE WEAK 
COUPLING LIMIT: THE SUBLEADING 
APPROXIMATION 

In the present section we demonstrate that in the sub- 
leading weak binding approximation, e — > 0, the parity 
degeneracy of the ground state obtained in the previous 
section is lifted: the parity splitting scales as A cx e 2 . Be- 
cause of this scaling the splitting is very small compared 
to the binding energy, A<6. 

Signs of the hole wave functions ip a ^ and 4>b^ given by 
Eqs. (|16p arc shown in Fig. [4j In this case the hole is 
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FIG. 4: (Color online). Left: signs of the wave function ip a -[, 
Right: signs of the wave function tp^ . The potential center is 
shown by the filled red circle. 

moving on the 'up' sublattice. To avoid misunderstand- 
ing we stress that there is only one hole, we do not show 
spins up in Fig. [Hjust to make the figure less busy. The 
figure clearly demonstrates that the states have opposite 
parities and different diagonal momenta, q a = (02, 02) 
and qi, = (02, —n/2). 

It is clear from Fig. |4] that the difference in energy 
between states ip a -[ and ipbt arises due to diagonal hop- 
ping of the hole in the vicinity of the potential, A cx 
t' e j j:\x(0)\ 2 , where t' e ^ is the effective diagonal hopping 
that is due to the bare if and also due to higher orders 
in t, see e.g. Eq. ([?])■ Moreover, the splitting cannot be 
just proportional to |x(0)| 2 ; the splitting must contain a 
gradient of \ because there is no a splitting for free hole 
propagation when \ = const. The first power of the gra- 
dient in the energy splitting is forbidden by parity. Thus 
we come to the following formula for the energy splitting 
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(19) 



(18) 



The formula contains the second power of gradient, so it 
is allowed by parity. Having in mind Eqs. p3|) . Q14p and 
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using Eq. (fl9|) wc conclude that 

A a tV«V^. (20) 

We stress that this formula follows from general symme- 
try considerations based on degeneracy of the free hole 
dispersion at the four nodal points q = (±7r/2, ±7r/2). 
The symmetry arguments certainly do not allow to de- 
termine a coefficient in Eq. (|2"D|) . However, they do allow 
us to determine the scaling law given by l[2"U)l . 

It is helpful to support the general considerations pre- 
sented in the previous paragraph by a numerical calcu- 
lation. Such a calculation in the regime t > J is hardly 
possible. However, in the regime t < J the calculation 
can be performed using results of Ref summarized in 
Eq.flT]). According to the results the spin quantum fluc- 
tuations can be integrated out and the hole propagation 
on the sublattice up is described by the following effective 
Hamiltonian 

H *ff =t'effJ2 h \ h i' + l "f f E h i h J" ■ ( 21 ) 
(ij>) (%]") 

Here h\ is the holon creation operator on the site i; all 
the sites i, f and j" belong to the sublattice up. To 
be specific we consider here the Coulomb attraction (jHJ), 
(|T0|> . The attractive interaction written in terms of holon 
operators reads 

He = -J2 Uih i hi ' ( 22 ) 

i 

where f7j is given by Eg. (fit))) . The Hamiltonian i? e // + 
He can be easily diagonalizcd numerically on a very large 
cluster. Results of diagonalizations for 30x30 cluster 
with t' e fj = 0.1, t'lff = 0.25 and for three values of the 
dimensionless charge Q = 0.75, Q = 0.5, and Q = 0.25 
are shown in Fig. O In this figure we show the holon 
probability distribution for shallow bound states, and in 
the legends wc present values of the rms radius, r rms \ the 
binding energy, e; and the ratio A/e 2 . According to the 
data in Fig. in the limit e — > the parity splitting A is 
decaying even slightly faster than oc e 2 . Most likely the 
small deviation from the expected e 2 scaling is due to the 
finite cluster size. We have also checked that the split- 
ting A vanishes at t' e ^ = 0. Altogether the numerical 
results presented in Fig. [5] confirm the scaling law given 
by Eq.QJJJ. 

The conclusion of the present section is that the ground 
state parity splitting is decaying A oc e 2 oc l/?y ms when 
the binding energy is decreasing, e — > 0, r rms — > oo. To 
estimate the coefficient in this dependence at t > J one 
can refer to results of exact numerical diagonalizations 
presented in Table I. It is known experimentally that in 
very lightly doped L^-zSr^CuC^ a hole binding energy 
to Sr ion is about e w — 10 meV, the bound state 'wave 
vector' is 0.4, and the rms radius of the bound state 
is r rms w 3, as discussed in Ref£ Estimates based on 




FIG. 5: (Color online). The shallow Coulomb bound state 
wave function squared, \ 2 > versus radius. Wave functions 
are shown for for three values of the dimensionless charge: 
Q = 0.75, Q = 0.5, and Q = 0.25. Values of the effective 
hopping parameters are t'^ii =0.1, t"tt = 0.25. In the legend, 
for every value of Q we also present the rms radius of the 
bound state, r rma ; the binding energy, e; and the ratio of the 
parity doublet splitting over the binding energy squared, A/e 2 . 



results derived in Sections IV and VI show that the ex- 
pected parity splitting of the ground state in this case is 
a small fraction of 1 meV. Therefore, parity breaking is 
practically a zero mode of the system. 



VII. PARITY BREAKING AND FORMATION 
OF THE LOCAL SPIN SPIRAL 

According to the discussion in previous sections a sin- 
gle hole bound state in the t- J model always has a definite 
value of the spin projection on the direction of staggered 
magnetization, S z = ±1/2, and it always has a definite 
parity. Dependent on parameters, t, t', etc, the ground 
state parity can be positive or negative, but it is definite. 
There is no local spin spiral at this stage. Very close 
to the ground state there is always a state of opposite 
parity. Wave functions of these states are given by Eqs. 
([H]) and (fTT|) . and parities are given by Eq. (fT5|) . Now, 
following the 2s-2p hydrogen atom scenario discussed in 
Section II, we can mix the opposite parity states by a 
weak external perturbation. There are two possibilities: 
l)mixing of states with different S z that belong to the 
same hole pocket, 2)mixing of states with the same S z 
that belong to different hole pockets. In the present sec- 
tion we consider the first possibility that leads to forma- 
tion of a local spin spiral shown in Fig. [1] The second 
possibility could lead to formation of a CDW. However, 
we show in the following section that this possibility is 
energetically unfavorable. 

Thus, let us mix the opposite parity states with dif- 
ferent S z that belong to the same pocket. To do so, we 
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impose a very weak uniform spin twist on the system. At 
this stage it becomes convenient to use the notation of 
the non-linear cr-model. In this notation the unit vector 
ft(r) shows direction of staggered spins. In the antifer- 
romagnetic state the spins are directed along the z-axis 
in the spin space, ft = n z = (0,0, 1). The uniform spin 
twist means that the spin direction ft rotates around a 
unit vector £ that is orthogonal to the z-axis. So locally 
we can write 



5ft(r) = (Q ■ r)[£ x ft] 



(23) 



Here Q -C 1 is the wave vector of the imposed twist. 
Let us direct Q along the b nodal direction, see Fig. [21 
Q = Qe^, where e& is the b- nodal unit vector. It is worth 
noting that generally directions in spin space and direc- 
tions in the coordinate space are completely independent . 
The interaction of a hole with the deformation of the spin 
fabric is of the following form^ 



Hi, 



-\piga x [ft X (e ■ Vn] 



(24) 



where e is a nodal unit vector corresponding to the par- 
ticular hole, a is the Pauli matrix acting on pseudospin 
of the hole. Note that in the notation of the original t- 
J model the effective Hamiltonian (|2"4")) is just the usual 
hole-spin- wave vertex shown in Fig.[SJ Therefore the cou- 




FIG. 6: Hole- spin-wave vertex corresponding to the Hamilto- 



pling constant is g « Zt « 1, where Z w J/t is the 
quasiparticle residue of the hole. 

Since Q = Qeb, the interaction (|2~4")) docs not mix 
states ipai and tp a i however, it does mix states ipbi and 
V>bj.- The corresponding interaction energy is 



5E int = -\/2<7Q(^|(f-a|V>> 



(25) 



If the interaction energy is larger than the parity doublet 
splitting, 



V2gQ > A , 



(26) 



the bound state wave function becomes a mixture of the 
opposite parity states 



(27) 



with the phase a determined by the condition (ip\a\ip) = 
£. Thus, the uniform spin twist Q is completely anal- 
ogous to a weak uniform electric filed E ext applied to 



hydrogen atom as has been discussed in Section II. The 
wave function mixing (|27[) is analogous to the mixing |T]) . 
Estimates based on values of A obtained in previous sec- 
tions show that for a bound state with radius r rms = 3 
the twist Q = 0.001 — 0.002 is already sufficient to break 
the parity according to Eqs. (f2"6")) . (|2~T|) . Note that this 
small value of Q corresponds to a wavelength of about 
5000 lattice spacing. 

The state (f27|) possesses a spin-flip dipolc moment and 
hence it creates a long range distortion of the spin fabric 
as has been discussed in Ref£ 



K x ft] 



9 (e 



[1 



e- 2Kr (l + 2 K r)] 



(28) 

Here p s « 0.18 J is the spin stiffness of the Heisenberg 
model, and k is the inverse radius of the charge core, see 
Eq. (fT3|) . Eq. (|28|) describes the local spiral depicted in 
Figd the local spin spiral is fully analogous to the long 
range scalar potential fi n d(f) generated by an excited 
hydrogen atom in a tiny external electric field, see Eq. 

To derive Eq. ([2"T]) and hence to justify the local spin 
spiral (|28p we have introduced a tiny external spin twist 
that enforces the parity breaking. Alternatively, one can 
consider an interaction between two holes bound to two 
impurities separated by a large distance r. Then there 
is no need for any external twist. Spin spirals induced 
by different holes lock each other. Hence the spin spiral 
induced hole-hole interaction isii 



g 1 1 

Eg ~ - 

4p s r z 



(29) 



This formula is valid at r < ta , while at r > r a the inter- 
action is Es oc 1/r 4 . Once more, this is absolutely similar 
to the case of two Hydrogen atoms, sec Eq.©. Estimates 
based on values of A obtained in Sections IV,VI show 
that for bound states with radius r rms = 3 the value of 
the crossover distance is ta ~ 50. So practically Eq. (|2"9"]) 
is always valid. To restore dimension in (j2"9")) one has to 
recall that g « J « 140meV, p s « 0.18J while dimcn- 
sionlcss r is expressed in units of lattice spacing. 



VIII. PARITY BREAKING AND POSSIBLE 
FORMATION OF THE CHARGE DENSITY 
WAVE 

We consider now a possible mixing of the bound states 
with the same S z that belong to different hole pockets. 
Since the spin is not changed there is no deformation of 
the spin fabric, and a usual electrostatic potential can 
mix the states. However, the spatial wave functions from 
different pockets differ by momentum K = (7r,0), or 
K = (0,7r). Therefore, to generate the mixing the elec- 
trostatic potential must be modulated at this momen- 
tum. So, the mechanism can produce a CDW with the 
wave vector K . Let is a Fourier component of the ex- 
ternal electrostatic potential. The component interacts 
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with the corresponding matrix element of charge density 
Pk = J ^ T (r)e- ife -> aT d 2 r = J X 2 (r)e^ K -^ r d 2 r 



[4k 2 + (K - fc) 2 ] 3 / 2 ' v ' 

We have used here Eqs. (fT6|) and (fl3|) . We assume that 
Pi =02, this allows to evaluate the integral in (|30|) an- 
alytically. Numerical integration shows that Eq. (|3"0")) is 
approximately valid even with non-equal inverse masses. 
For example at f3i //?2 = 4 the deviation from the analyt- 
ical expression Eq. (f3"0")) does not exceed a few per cent. 

We proceed now directly to the Coulomb interac- 
tion between two holes bound to two different impuri- 
ties separated by large distance r. Since the system is 
two-dimensional the electrostatic potential created by a 
charge density component pk is 

2tt . 

4>k = -jrPk ■ (31) 

Therefore the Coulomb interaction energy between two 
spatially separated bound states reads 

*>-!* (32 » 

Note that the sign is negative because the system always 
tunes up the mixing phases to reduce energy. Evaluation 
of ([32)) with account of (|30|) is straightforward. In the 
limit Kr > 1 the result reads 

* - - © ■ 

Here n is dimcnsionless and we put the factor e 2 / (ae e ) m 
95meV to restore the dimension of energy (e is the el- 
ementary charge, a = 3.81 A is the lattice spacing and 
e e w 40 is the dielectric constant). 

Now we can compare the CDW Coulomb interaction 
(|33|) with the spin-spiral interaction ([29]) . In LSCO the 
'wave vector' of the bound state is k w 0.4, see Ref.— 



Let us take r = 4 that corresponds to the average dis- 
tance between bound states at the doping level x f» 0.06. 
With these parameters one finds Ec < lmeV while 
Es ~ 15meV. Thus formation of the CDW is energet- 
ically unfavorable compared to formation of the spin spi- 
ral . 



IX. CONCLUSIONS 

We have considered a single hole injected into a two 
dimensional Mott insulator on a square lattice with a 
long range antiferromagnctic order. The system is de- 
scribed by the extended t-J model. An important point 
is that minima of the hole dispersion are at nodal points 
(±7r/2, ±/tt/2). The hole is bound by an impurity po- 
tential. The impurity is located at a center of the lattice 
plaquette, so the potential itself does not break the local 
square lattice symmetry. 

1) All bound states have definite parity and they are 
doubly degenerate with respect to the spin projection on 
the axis of the staggered magnetization, S z = ±i. 

2) The ground state always has a very close state of op- 
posite parity (parity doublet). For shallow bound states 
splitting within the parity doublet scales as A cx e 2 , 
where e is binding energy. 

3) For shallow bound states the parity splitting A is 
extremely small. Therefore an extremely small external 
twist of the spin fabric breaks parity. The breaking cre- 
ates a long range spiral distortion of the spin fabric. The 
breaking can be also created by another impurity; in this 
case the local spirals of two impurities lock each other. 

4) The bound state parity breaking in the t-J model 
is very similar to the parity breaking within the 2sx/2 — 
2pi/2 parity doublet of the hydrogen atom. 
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